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The task of decoupling, i.e., removing unwanted interactions in a system Hamiltonian and/or 
couplings with an environment (decoherence) , plays an important role in controlling quantum sys- 
tems. There are many efficient decoupling schemes based on combinatorial concepts like orthogonal 
arrays, difference schemes and Hadamard matrices. So far these (combinatorial) decoupling schemes 
have relied on the ability to effect sequences of instantaneous, arbitrarily strong control Hamiltoni- 
ans (bang-bang controls). To overcome the shortcomings of bang-bang control Viola AND Knill 
proposed a method called "Eulerian decoupling" that allows the use of bounded-strength controls 
for decoupling. However, their method was not directly designed to take advantage of the composite 
structure of multipartite quantum systems. In this paper we define a combinatorial structure called 
an Eulerian orthogonal array. It merges the desirable properties of orthogonal arrays and Eulerian 
cycles in Cayley graphs (that are the basis of Eulerian decoupling). We show that this structure 
gives rise to decoupling schemes with bounded-strength control Hamiltonians that can be applied 
to composite quantum systems with few body Hamiltonians and special couplings with the environ- 
ment. Furthermore, we show how to construct Eulerian orthogonal arrays having good parameters 
in order to obtain efficient decoupling schemes. 

PACS numbers: 03.67.Lx, 03.65.Fd, 03.67.-a 



I. INTRODUCTION 

An important task in quantum information theory 
consists in selectively removing unwanted contributions 
of the system Hamiltonian and/or switching off cou- 
plings of the quantum system to an uncontrollable 
environment (the later being responsible for decoher- 
ence). This ta s k is usually cal led decou pling (see e.g . 
JvLM IVKLQfl IVLK991 IZanQiH and e.g. |JKM ISMOl 
iLCYYOft IWR..IB02L iLeuO^ for schemes using combi- 
natorial concepts). More generally, one is also inter- 
ested in effectively changing the system Hamiltonian 
in order to simulate some desired Hamiltonian; this is 
usually referred to as simulating Hamiltonians (see e.g. 
jW.TBnl IDNBT021 IWR.TBOl lBCL+02| ). In this pa- 
per we will concentrate on designing efficient decoupling 
schemes. 

Methods of dynamical decoupling and also simulat- 
ing Hamiltonians derive their basic physical intuition 
from coherent averaging techniques in high-reso lution nu- 
clear ma gnetic resonance (NMR) spectroscopy |WHH68. 
EBW87J. A decoupling scheme is understood as a con- 
trol protocol which relies on the repeated application of 
controls drawn from a finite set in order to change effec- 
tively the natural time evolution to the desired evolution. 
Many efficient decoupling schemes can be designed with 
the help of combinatorial concepts like e.g. Hadamard 
matrices, difference schemes and orthogonal arrays. The 
entries of these structures describe how to choose the con- 
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trols. The reason why it is possible to use these combina- 
torial objects is the special structure of the system Hamil- 
tonians (pair-interactions or more generally few body 
Hamiltonians). So far all these (combinatorial) schemes 
relied on the ability to effect sequences of instantaneous, 
arbitrarily strong controls Hamiltonians (bang-bang con- 
trols). Because of the requirement of bang-bang controls 
such schemes are unrealistic in many situations. To over- 
come the shortcomings of bang-bang decoupling Viola 
AND Knill proposed a general method (called Eulerian 
decoupling) for implementing decoupling with bounded 
controls (i.e., continuo usly modulate d bounded-strength 
control Hamiltonians) |VK02llVio04j . This method offers 
many advantageous over bang-bang decoupling. How- 
ever, their method was not directly designed to make use 
of the special structure of few body Hamiltonians to re- 
duce the complexity of decoupling. 

We show how to incorporate some of the above com- 
binatorial methods (that were used so far only in the 
bang-bang formulation) into the method by Viola and 
Knill in order to obtain efficient decoupling schemes 
with bounded controls. Our schemes can be applied to 
few body Hamiltonians and special couplings with an 
environment. Our schemes rely on a combinatorial ob- 
ject (which we call) an Eulerian orthogonal array. We 
show how to construct these objects with good parame- 
ters from error correcting codes. 

The paper is organized as follows. In section 2 we 
describe the principles of dynamical decoupling. We de- 
scribe briefly the so-called first-order approximation that 
is the basis for all decoupling schemes. In section 3 we 
recall how to construct decoupling schemes with bang- 
bang and bounded strength controls. The method using 
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bounded strength control is Viola AND Knill's Eule- 
rian decoupling. Here no special structure of the quan- 
tum system is assumed. Then in section 4 we consider 
quantum systems consisting of coupled qudits. We first 
recall in subsection 4.1 how to construct efficient decou- 
pling schemes with bang-bang controls with the help of 
orthogonal arrays. Efficiency means that the number of 
necessary pulses grows polynomially with the number of 
qudits. In subsection 4.2 we show how to merge the con- 
cept of orthogonal arrays with the idea of Eulerian de- 
coupling in order to obtain efficient decoupling schemes 
with bounded controls. Our method is based on a com- 
binatorial structure called Eulerian orthogonal array. In 
sections 5 we show how to construct Eulerian orthogonal 
arrays with good parameters. 



II. PRINCIPLES OF DYNAMICAL 
DECOUPLING 

A decoupling scheme is understood as a control pro- 
tocol which relies on the repeated application of controls 
drawn from a finite set in order to change effectively the 
natural time e volution to th e desired evolution. We refer 
the reader to |VK02llVio04| for a more detailed descrip- 
tion. In the following we give a brief introduction based 
on the above articles. 

The joint evolution of the target system S in interac- 
tion with the environment E is described by a total drift 
Hamiltonian of the form 



H = H s ® 1 E + 1.' 



Hp 



Use , Hse — Sa® E a , 



. (1) 

where Hs and He characterize the isolated dynamics of 
the system and the environment, respectively, and the 
interaction term Hse is responsible for introducing un- 
wanted decoherence effects and dissipation effects in the 
reduced dynamics of S alone. Without loss of general- 
ity we will always choose the operators S a and Hs to be 
traceless. 

The idea behind dynamical decoupling is to add a spe- 
cially designed controller, described by a time-dependent 
control Hamiltonian H c (t) acting on only the target sys- 
tem S 1 , in such a way that the resulting controlled dynam- 
ics is described by an effective Hamiltonian H c g which no 
longer contains any coupling terms between S and E, i.e., 



H eS =H s ®l E + ls®H E , 



(2) 



for an appropriate, possibly modified, system Hamilto- 
nian Hg. In this paper we will be interested in the case 
that H s = 0. 

Decoupling protocols are most conveniently con- 
structed by directly looking at the control propagator 
associated to H c (t) 



U c (t) =Texp|-i^ ^ c (r)rfr| , 



(3) 



where T denotes the time ordering. 

Decoupling is based on the so-called first order de- 
coupling. The control actions are always cyclic, i.e., 
U c (t + T c ) = U c (t) for some cycle time T c and for all 
t. The stroboscopic dynamics U{pia) with tu = MT C 
and M S N may be described by a propagator 



U(t M ) = exp(-iHt 



M) 



(4) 



for a time-independent effective Hamiltonian iJ c g . If, in 
addition, T c is sufficiently short, then the effective Hamil- 
tonian is accurately represented by the following lowest- 
order Hamiltonian 



H 



1 



Ul{t)HU c {t)dr . 



(5) 



While higher-order terms can be systematically evalu- 
ated, the approximation in JSJ) becomes more and more 
exact as the fast control limit T c — > is approached. 
First-order decoupling is based on this approximation. 

Having introduced the framework of decoupling we 
address the problem of designing efficient decoupling 
schemes with first bang-bang and then bounded-strength 
controls for general Hamiltonians and couplings with the 
environment. 



III. DECOUPLING SCHEMES 

In this section do not assume any special structure 
of the target system S, the Hamiltonian Hs, and the 
coupling to the environment given by S^s. We first 
discuss how to realize decoupling with controls of un- 
bounded strength (bang-bang) and then with bounded 
con trol (E ulerian decoupling) . The presentation is based 
on |VK02j . 



A. Bang-bang control 

The time-average in JSJ can be expressed directly as an 
average over a group in the following simple bang-bang 
decoupling setting. Let G be a discrete group of order 
\G\ acting on the Hilbert space of the target system Tis 
via a faithful, unitary, projective representation, 



G 
9 



U{H S ) 

P(g) 



u„ 



(6) 



where U(Hs) denotes the group of unitary matrices act- 
ing on Hs- Let A G N and A > 0. 

A decoupling scheme using U g 's as control operations 
is specified by a sequence (51,32, ■ ■ ■ >9n) with N := A|G| 
and entries from G. The number N is called the length 
of the decoupling scheme. The entries g/s specify the 
control propagator U c (t) over each of the A|G| equally 
long subintervals. A control cycle is defined by 



U c ((j-l)A + r) = U gj , re[0,A) 



(7) 
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with T c = \\G\A for some A > 0, and j = 1, . . . ,N. 

If all group elements appear exactly A times in the list 
(<7i, . . . , pjv) then the resulting control action corresponds 
to extracting the G-invariant component of X. We have 

1 N A 

- £ ul jXUgj = - £ ulxu g = U G (X) , 



where 



j=i sec 



1 1 5£G 



(8) 



Note that if the representation p in © is irreducible than 
we have Hc(X) = tr(X)/d Id for all X, where d is the 
dimension of Us- 

An example for such an irreducible, unitary, projec- 
tive representation is given in the following. The discrete 
Fourier transform of length d G N is the unitary transfor- 
mation defined by DFT^ := Y^Te^o^ '*l^)(^l> wnere 
u> denotes the primitive d-th root of unity e 27ri / d . Next, 
define operators S := Y^t=o + M- where the indices 
are reduced modulo d, and T :— DFTt • S ■ DFT d = 



J2 k Z v k \k)(k\. Then the map 



P 



Z d x Z d 
( 0> 6) 



W(d) 



(9) 



is an irreducible, unitary, projective representation. Note 
that for d = 2 one obtains 1 , a x , a y , cr z , where the cr's are 
the Pauli matrices. 

Now it clear that we can remove the couplings with the 
environment and switch off the natural time evolution 
of the quantum system by performing the control opera- 
tions according to and the representation in ©. This 
is because Hg(So,) = for all a and Hg(-Hs) = 0. How- 
ever, this method has the following disadvantage that 
makes it unrealistic form in many situations. Accord- 
ing to the rule in Q the control propagator U c {t) jumps 
from U 9j _ 1 to U gj = (U gj U g ~._ 1 )U gj through the applica- 
tion of an arbitrarily strong, instantaneous kick at the 
jth endpoint tj = jA, realizing the bang-bang pulse 
U Sj = Ug j U^._ 1 with Sj := g^gTi (equality is under- 
stood here up to a phase factor). In the next section we 
describe how to avoid such bang-bang controls. 



B. Eulerian decoupling 

As already mentioned in the introduction the require- 
ments for bang-bang control are highly unrealistic. Vi- 
ola and Knill proposed a method called Eulerian de- 
coupling |VK02l ] that avoids the use of such bang-bang 
pulses. In Eulerian decoupling the control propagator 
U c (t) is varied smoothly from U gj _ 1 to U 9j by a control 
action distributed along the whole jth subinterval. 

Let 5* be a generating set for G, i.e., any element of 
G can be written as a product of elements of S. The 



Cayley graph T(G, S) of G with respect to S is a directed 
graph whose vertices are labeled by the group elements 
and whose edges are labeled by the generators. More 
precisely, the vertex g is joined to the vertex h if and 
only if gbr 1 — s for some s € S, i.e., g = sh. 

It is assumed that we have the ability to physically 
implement the generators s £ S, i.e., to implement the 
unitaries U s by the application of some suitably chosen 
control Hamiltonians h s (t) over A: 



U s =u s (A) 



where 



u s (5) 



T | exp i 



,{r)dr) 



(10) 



(11) 



for (5 G [0, A]. The choice of the control Hamiltonians 
h s (t) is not unique. This allows for additional flexibility 
for the concrete implementation. Once a choice of the 
control Hamiltonians is made, the control action is deter- 
mined by assigning a cycle time and a rule for switching 
the Hamiltonians h s ( t) durin g the cycle subintervals. 

Viola and Knill |VK02| showed that decoupling can 
be achieved by sequentially implementing generators so 
that they follow a Eulerian cycle in T(G, S). An Eule- 
rian cycle is defined as a cycle that uses each edge exactly 
once. Because a Cayley graph is regular, it always has an 
Eulerian cycle, whose le ngth is necessarily V = \G\ \S\ 
(see e.g. |Bol98t lGR0ll | for the definition of these no- 
tions). For our purposes, we use a slightly more general 
definition: an Eulerian cycle with multiplicity A is a cy- 
cle that uses each edge exactly A times. Clearly, such 
an Eulerian cycle has necessarily length V = A \ G\ \S\. 
We will choose an Eulerian cycle to begin at the identity 
element of G. Therefore, an Eulerian cycle can be de- 
scribed as a list (si, . . . , sn) with entries from S. Each 
entry identifies the edge via which we leave the vertex. 

Decoupling according to an Eulerian cycle C :— 
(si, . . . , sn) is defined by setting the cycle time T c = NA 
and by choosing the control propagators U c (t) as follows: 



U c ((j - 1)A + S) = u Sj (S) U c ((j - 1)A) 



(12) 



where S e [0, A) and u Sj (6) is defined in (fTOf) and ltTT|l . 
This decoupling prescription means that during the jth 
subinterval one choses as a control Hamiltonian the one 
that realizes the generator Sj, i.e., the jth element of C. 

The effective Hamiltonian H under Eulerian decou- 
pling is obtained by evaluating the time-average in (jSJ) 
with the control propagator being given by (|12fl . The 
resulting V terms can be partitioned in \S\ families, each 
corresponding to a fixed generator. Because for each s 
the cycle C contains exactly A s-labeled edges ending at 
any given vertex g, each family leads to a sum over the 
group elements as in (JSJ. 

For these reasons the quantum operation Qc defined 
by C can decomposed as 



Q C (X) =U G (F S (X)) 



(13) 
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with the map F$ implementing an average over both the 
group generators and control sub-interval: 

Fs(X) = |1 £ Mr)H S )Xu s (r)dr . (14) 

The link between Eulerian decoupling and bang-bang de- 
coupling by averaging over G is established in the follow- 
ing theorem. Some additio nal compatibility between ILq 
and Fs is necessary |VKf)2j . Let us repeat all the notions 
before stating the theorem. Let G be a group that acts 
via a faithful, unitary, projective representation g i— > U g 
on C d . The decoupling group algebra V of G is the C- 
linear span of the matrices U g . 

Theorem 1 (Eulerian decoupling) 

Let X be any operator acting on C d . If the control Hamil- 
tonians h s (t) are in the decoupling group algebra, i.e., 
h s (S) £ V for all times S £ [0, A] and all s £ S, then Eu- 
lerian decoupling according an Eulerian cycle C as spec- 
ified by the the rule in f Wjl has the same effect as aver- 
aging over G as in i.e, 

Qc(X) = U G (X). 

For the proof we refer the reader to VK02J. Note that the 
bang-bang limit is formally recovered by substituting the 
map Fs by the identity map. In the Eulerian approach, 
at the expense of lengthening the control cycle, the same 
G-symmctrization can be attained using only bounded- 
strength controls. The maximum strengths achievable 
in implementing the generators directly bounds the min- 
imum attainable T c , and therefore the accuracy of the 
first-order approximation. 

IV. EFFICIENT DECOUPLING SYSTEMS 

In this section we consider a target systems that is 
composed of n coupled qudits, i.e., its Hilbert space TLs 
is given by the tensor product TLs — (C d ) <8 ™ . We say 
that a family of decoupling schemes is efficient if the 
number of control operations grows polynomially with 
the number of qudits. So far there were only efficient de- 
coupling schemes using bang-bang controls (see e.g. the 
references given in the introduction). The schemes rely 
on the special structure of the system Hamiltonians and 
the couplings to the environment. It is assumed that 
the system Hamiltonian is a so-called few body Hamilto- 
nian. To define this precisely we need to introduce some 
notions. For any operator A acting on C d we denote by 
the operator that acts as A on the kth qudit, i.e., 
A^ = l<g)"-®l(g>A(g>l<g)"-<g>l. 

Let B :— {a a \ a = 1, . . . , d 2 } be an basis of for the 
vector space C dxd of matrices acting on <C d . We say 
that an operator X acts on the qudits k\ , . . . , k t with 
1 < k\ < . . . < fc{ < n if it can be expressed as follows 

X- V r a (k ^ ■ ■ ■ cr {kt) 

oil,..., at 



for some x aii ... t0tt S C. 

We assume that the system Hamiltonian is a i-body 
Hamiltonian, i.e., it can be decomposed as 

H S := J2 H *,:.te> (15) 

ki,...,kt 

where Hk t ,...,kt are traceless operators acting on qudits 
k\ , . . . , kt only. For t — 2 one also says that H is a 
pair-interaction Hamiltonian. Furthermore, we assume 
that the environment couples independently to t-tuples 
of qudits, i.e., we have 

Hse= Yj S ki,...M®E ku ... M (16) 

k\ ,...,kt 

where Sk lt ...,k t are traceless operators acting on qubits 
fci, . . . , kt only and i?/ Cl ,...,fc t act on the Hilbert space He 
of the environment. 

It will be convenient to use the following definition. Let 
X be an arbitrary operator acting on (C d )®*. We define 
its embedding x'* 1 '— •**) into (C d )® n to be the operator 

X (ku.-.M) = J- x ai ,... at a^---a^\ (17) 

a\,...,ott 

where X = a( x ail ,„ at a ai ® ■ • ■ l S>cr at is the expan- 
sion of X in the product basis B® 1 . 

A. Decoupling with bang-bang controls based on 
orthogonal arrays 

We assume that we can perform bang-bang controls on 
each qudit individually. Formally, all control operations 
are elements of some finite subset of the group U{cS)® n , 
where U(d) denotes the group of unitary matrices acting 
on C d . In the following we recall how orthogonal arrays 
may be used to construct efficient decoupling schemes. 
Orthogonal arrays appeared first in statistics where they 
were used in the design of experiments for collecting sta- 
tistical d^^wstemAticaJh^. We refer the reader to the 
books |BJL99t ICD9l lHSS99j for applications and con- 
structions of orthogonal arrays. Stollsteimer and Mahler 
first used orthogonal arrays (or OAs for short) for the 
construc tion of decoupling schemes and selective coupling 
schemes ;S.\1() \\ for qubit systems with pair-interactions. 
This method was g eneralized to qudi t systems with t- 
local interactions in |WRJB02l lRW04j . 

Definition 1 (Orthogonal array of strength t) Let 

A be a finite alphabet and let n, N G N. An n x N array 
M with entries from A is an orthogonal array with \A\ 
levels, strength t, and multiplicity A if and only if every 
t x N sub-array of M contains each possible t-tuple of 
elements in A 1 precisely A times as a column. We use 
the notation OA\(N,n,s,t) to denote a corresponding 
orthogonal array. If X, s, and t are understood we also 
use the shorthand notation OA(N, n) . 
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An important special case arises if the strength t is two. 
This means that each pair of elements of A occurs A times 
in the list ((akj,aij) j — 1, . . . N) for 1 < k < I < n. 
Most of the known c onstruction actually yield arrays of 
strength two HSS99]. For many physical systems it will 
be sufficient to study arrays of small strength since the 
strength relates to the degree of the interactions, i.e., 
for pair-interaction Hamiltonians it is sufficient to con- 
sider arrays of streng th t = 2. For an example of such 
orthogonal arrays see |Rot04lRWf)4| . 

The basic idea is to use an orthogonal array M with 
parameters OA(N, n, d 2 , 2) over an alphabet A of size 
d 2 . Here d denotes the dimension of the qudits. The ele- 
ments of A are identified with the elements of the group 
G := Zd x Zd that acts irreducibly on C d via the map in 
©. The columns (gxj, . . . ,g n j) T of M specify the con- 
trol propagators U c (t) over each of the N equally long 
subintervals. A control cycle is defined by 

U c ((j - 1)A + r) = U m ®U g2] ®---® U 9nj , (18) 

where r e [0, A), T c = NA for some A > 0, and j = 
1,...,N. 



The first equality is because X^ kl, '"' kt ' acts on the qudits 
ki, . . . ,k t only. Note that the representation of G xt to 
11(d)®* given by (hi, ...,h t ) U hl ® • ■ ■ ® U ht is irre- 
ducible. Since M is an orthogonal array of strength t the 
list ((gki,j, ■ ■ ■ >9k t ,j))f=i contains every element of G xi 
exactly A times. Therefore, we average over the group 
G xt acting via an irreducible representation as in JSJ. 
This proves eq. 1(277)1. □ 



It is clear that we can switch off all Hamiltonians of 
the form in (|15fl and all couplings with an environment 
of the form in lltj|) with control actions specified by the 
above theorem. This is because all operators Hk lt ...,kt 
and Sfa ...,fc t have the form X ( fe i. ■■■■'=*) for some traceless 
operator X acting on (C 1 )®*. 



The following theorem shows that the prescription in 
(11811 allows to decouple few body Hamiltonians and cou- 
plings with the environment. Let G be an arbitrary finite 
group. We denote by G xt the direct product G x • • • x G 
having t components. 

Theorem 2 (Decoupling with OAs) 

Let G := Zd x Zd and g i— > U g be the irreducible, unitary, 
projective representation as in Jty. Let M = (gkj) be an 
OA(n, N) be an orthogonal array of strength t over the 
group G. Let Um denote the control action that corre- 
sponds to $181) . Then we have 

n M (A (fcl '-< fct) ) = (19) 

for an arbitrary traceless operator acting on (C d )®* and 
any t-tuple (hi, ■ ■ ■ , kt) with 1 < k\ < . . . < k t < n. 

Proof. The idea is to reduce the problem for each t- 
tuple (ki, . . . ,k t ) to the case in ® by using the special 
structure of the operator X^ kl ''"' kt ' . We have 




-i (fci,...,fc t ) 



- (fei,-,fe t ) 

• ® u ht ) 

(20) 



B. Decoupling with bounded strength controls 
based on Eulerian orthogonal arrays 

Finally, we show how to combine the ideas of Eule- 
rian decoupling and orthogonal arrays. This is done by 
introducing the concept of Eulerian orthogonal arrays. 

Definition 2 (Eulerian orthogonal array) 

A n x N -matrix M — (gkj) with entries from the group G 
is said to be an Eulerian orthogonal array of strength t iff 
for all t-tuples (k\, k^,. . . ,kt) with 1 < ki < . . . < kt < n 
there is a generating set Sk lt k2,...,kt of G such that the 
list of group elements 

{(9kuj,9k 2 ,j, ■ ■ ••'//. .,) ( V , (21) 

defines an Eulerian cycle in the Cayley graph 
T(G xt , 5fe 1 ,fe 2 ,...,fe t ). 



n M (* (fel --' fct) ) 



N 



i N 



■u g „y x 



u Bhn ^ x{u ghi .®u gk2 .® 



J2 ® U h2 ® • • • <8> U ht ) ] X (U ht ® U h2 

(h 1 ,...j lt )eG>< t 

[n G x t (A)] (fel — kt) = o( fci — fct ^ = 0. 

I 
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Note that the above conditions automatically implies 
that M is a (usual) orthogonal array of strength t. 

We assume that we have the ability to implement the 
group elements g G G, i.e., to implement the unitaries 
U g on the individual qudits by the application of control 
Hamiltonians h g (t) over A as in Hl()|) and JUJ. This 
means that we have the ability to switch on the control 
Hamiltonians h g (t) on any qudit, i.e., 1(g)- • -® l®fo ff (i)® 
1 ® •••(8)1. 

We define decoupling according to an Eulerian orthog- 
onal array M = (gkj) by setting the cycle time T c = NA 
and by assigning the control propagators as follows: 

U c ((j-l)A + 6) 
= (u Slj {5)®---®u Snj (5))U c ({j-l)A) (22) 

where S € [0, A) and s kj = g^gkj+i for j = 1, . . . , iV- 1 

and s k N = 9k N gk\- The tuples (sk t ,j, ■ ■ ■ , Sk t ,j) are the 
edges in the Eulerian cycle defined by the rows ki,...,kt 
of M. 

Theorem 3 (Decoupling with Eulerian OAs) 

Let G := x and g i— > U g be the irreducible, unitary, 
projective representation in Ity. Let M = (gkj) be an Eu- 
lerian orthogonal array over G of size nx N and strength 
t. Let Qm denote the control action that results from the 
control propagator defined in \2ty) . Then we have 

Q M (X^'-' k ^) = (23) 

for an arbitrary traceless operator X acting on (C d )® n 
and any t-tuple (k\, . . . , fc t ) with 1 < k\ < . . . < k n < n. 

Proof. Again the idea is to reduce the problem for each 
t-tuple (fci, . . . , kt) to the case of Theorem ^ by using 
the special structure of the operator X^ kl ''"' kt ' . Let us 
denote by Ck lt .... kt the Eulerian cycle in the Cayley graph 
r(G x *, Sk lr ..,kt) that is define by the rows ki,...,kt of 
M . Then we have 

Q M (X^-^) = [Q Cki fet ( X )] (fel - fct) (24) 

= [n GXt (F Sfei ,..., fct (A))] (fcl - fet) 

_ Q(fcl,...,fc t ) _ Q ^ 

Eq. |J2U is due to the fact that X^ 1 '-'^ acts on qu- 
dits fci, . . . , kt only. The remaining equalities follow from 
Theorem Q] because all its conditions are satisfied. □ 

Again it is clear that we can switch off all Hamiltonians 
of the form in l|15|) and all couplings with an environment 
of the form in i|16|) with control actions specified by the 
above theorem. 



V. EULERIAN OAS FROM LINEAR ERROR 
CORRECTING CODES 

We will construct Eulerian orthogonal arrays using lin- 
ear error correcting codes. Let us briefly repeat some 



basis facts about linear error correcting codes and their 
relationship to orthogonal arrays. 

A linear code over the finite field ¥ q is a fc-dimensional 
subspace of the vector space ¥ q . We consider finite 
fields of q = d 2 only. In this case the additive group 
of the finite field ¥ q is isomorphic to Z c i x Z c i\ we will 
again use the irreducible representation in J§J|. The 
space F^ is endowed with a metric called Hamming dis- 
tance. It is defined as follows: for x — (x%, . . . , x n ) G F™ 
we have that wt(x) := \{i G : Xi ^ 0}[. 

The minimum distance of a linear code C is defined by 
d — d min :— min{wt(c) : c G C, c ^ o}, where o denotes 
the zero vector. In this situation we say shortly that C 
is an [n, k, d] q code. We need the fact that a [n, k] q linear 
code can be described by a generator matrix G of size 
n x k with entries from ¥ q . The matrix G defines the 
embedding from F^ to F"; the code words c G C are the 
images of the vectors m G ¥ q , i.e., c = Gm. We need one 
more definition which is the dual code C 1 - of C defined 
by C x := {x G F™ : x ■ y = for all y G C}; the dot 
product x ■ y is given by ^Z i=1 x%yi. In the following we 
refer to the minimum distance d of the dual code as the 
dual distance. 

The following theorem HSS99, Theorem 4.6] estab- 
lishes a close relationship between orthogonal arrays and 
error-correctin g codes. This theorem was also used in 
|Rotn4lKW7)4| . 

Theorem 4 (OAs from linear codes) 

Let C be a linear [n, k, d] q code over¥ q with dual distance 
d . Arrange the codewords of C into the columns of a 
matrix A G ¥ q xqk . Then A is an OA(q k ,n,q,d 1 - - 1). 

Now we show how to modify the above construction in 
order to obtain Eulerian orthogonal arrays. 

Theorem 5 (Eulerian OAs from linear codes) 

Let C be a [n,k] q -code with dual distance d and G be 
a generator matrix for C. Let C := (mo, • • • , TOjv— i) 
be an Eulerian cycle in the Cayley graph T(V, S) with 
multiplicity 1, where the group is V := F^ and the 
generating set is the group itself, i.e., S := ¥ k . The 
length of such an Eulerian cycle is necessarily N = q 2k . 
Set t := d — 1. Then the n x N matrix M whose 
columns are defined to be Gnij for j = 0, . . . , N — 1 
is an Eulerian orthogonal array over ¥ q of strength t. 
Furthermore, we have Sk lt ...,k t = ^ Xt f or a ^ t-tuples 
(k%, . . . , k t ) with 1 < k\ < . . . < kt < n. 

Proof. Since C = (mo, . . . ,mjv-i) is an Eulerian cycle 
all elements of V = ¥ k appear exactly q k (corresponding 
to the size \S\ — q k of the generating set S) times in 
C. Therefore, the column vector Gm appears exactly q k 
times in M for all m €¥ k . It now follows from Theorem^] 
that M is an orthogonal array; its multiplicity is just q k 
times the multiplicity of an OA constructed based on 
Theorem 0J 
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Let to be an arbitrary element of V = F^ and I m := 
{j | 771 j = to}. Then every element of S = F^' appears 
exactly once in the list (rrij+i — raj \ j G I m ) because 
C is an Eulerian cycle in r(V, S) with multiplicity one 
(the addition is done modulo N). Consequently the list 
of transitions that occur in M from all columns of the 
form Gra, i.e., (Gmj+i — Grrij \ j G I m ) is independent 
of to and is equal (up to reordering the columns) to the 
orthogonal array M' :— (Ge | e G F£); it follows from 
Theorem |U that M' is an orthogonal array. This proves 
that M is Eulerian and also that 5fe 1 ,...,fc t = G xt for all 
t-tuples since M' is an orthogonal array of strength t. □ 

Note that our constr uction is closely related to 
Rotteler's construction |Rot04j . In that paper Hamil- 
tonian cycles in the Cayley graph r(Fg,5) are used, 
where the generating set is S given by the k coordinate 
vectors. The motivation behind this construction is to re- 
duce the number of different control pulses in bang-bang 
decoupling. 

Let us now explain how to construct decoupling 
schemes for general i-body Hamiltonians acting on n qu- 
dits with bounded controls based on Theorem |SJ To 
obtain a decoupling scheme using a minimal number of 
pulses we have to find a code [n, k] q such that k is mini- 
mal and the dual distance d 1 - is at least t + 1. This may 
be formulated in terms of the dual code which has pa- 
rameters C = [n,n — k, d^q. The dual code C ± should 
contain the maximally possible number of code words for 
given n and d 1 " . This question is one of the central opti- 
mization problems in the theory of error correcting codes. 
To find such optimal or best known codes one could e.g. 
use the computer algebra system MAGMA |BCP97t | that 
contains a table of best known linear codes (i.e., with 
the maximal number of code words) for given length and 
minimal distance. 

We now consider a quantum system consisting of n 
qubits which are governed by a pair-interaction Hamil- 
tonian. For such a system we can construct decou- 
pling schemes using N pulses from an orthogonal array 
OA(N, n,4, 2). Hence, in order to apply Theorem0]and 



[3] we have to find a code C of linear codes over F4 for 
which the parameters are [n, k, d] and for which the dual 
distance is at lea st 3. This can be done with the help of 
Hamming codes |Rot04 IRWO^ . For every to G N there 
is an orthogonal array OA(4 m ,(4 m - l)/3,4,2). The 
columns of this OA are codewords of the dual code of 
a Hamming code. The corresponding Eulerian orthogo- 
nal array has parameters OA(16 m , (4 m - l)/3, 4, 2). 

To obtain a decoupling scheme for a quantum system 
consisting of n qubits, where n is an arbitrary natural 
number, i. e., not necessarily of the form n = (4 m — l)/3 
we proceed as follows: first let to G N be the unique in- 
teger such that n < 4 ~ 1 < An. Then construct the or- 
thogonal array with parameters CL4(4 m , (4 m - l)/3, 4, 2) 
for bang-bang controls and the Eulerian orthogonal array 
with parameters OA(16 m , (4 m - l)/3,4,2) for bounded- 
strength controls, respectively. These results shows that 
the complexity of decoupling for general pair-interactions 
Hamiltonians acting on n qubits scales at most linearly 
in 71 for bang-bang controls and at most quadratically in 
n for bounded-strength controls, respectively. 

Conclusions and Discussions. — We have shown that 
it is possible to construct decoupling schemes using 
bounded-strength controls for composite multipartite qu- 
dit systems with the help of Eulerian orthogonal arrays. 
Our concept of Eulerian orthogonal arrays merges the de- 
sirable properties of usual orthogonal arrays (that were 
use to construct efficient decoupling schemes with bang- 
bang controls) and Eulerian cycles that are at the heart 
of Viola and Knill's Eulerian decoupling method. We 
have shown how to construct efficient Eulerian orthog- 
onal arrays based on linear error correcting codes. It 
would be interesting to find new construction of Eulerian 
orthogonal arrays that might yield decoupling schemes 
with a smaller number of pulses. 
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